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We present a theory for carrier induced ferromagnetism in diluted magnetic semi-conductor (DMS). 
Our approach treats on equal footing quantum fluctuations within the RPA approximation and dis- 
order within CPA. This method allows for the calculation of T c , magnetization and magnon spec- 
trum as a function of hole, impurity concentration and temperature. It is shown that, sufficiently 
close to T c , and within our decoupling scheme (Tyablicov type) the CPA for the itinerant electron 
gas reduces to the Virtual Crystal Approximation. This allows, in the low impurity concentration 
and low density of carriers to provide analytical expression for T c . For illustration, we consider the 
case of Ga\- c Mn c As and compare our results with available experimental data. 
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The discovery of carrier induced ferromagnetism in 
DMS have attracted considerable attention from both 
theoreticians and experimentalists. The interest for these 
material is mainly stimulated by the possible technolog- 
ical applications (e.g. semi-conductor spin devices). For 
example by doping GaAs [QJ^] with magnetic impurities 
Mn 2+ , T c exceeding 100 K has been reached. The dop- 
ing of a III-V semiconductor compound with Mn impu- 
rities introduces simultaneously local magnetic moments 
(S = 5/2) and itinerant valence band carriers (s = 1/2). 
One of the important open issues is to find out whether 
it is possible to reach critical Curie temperature of or- 
der 300 K. Thus it is important to understand theoreti- 
cally how T c varies with the impurity concentration, ef- 
fective mass, hole concentration and exchange integral. 
Many theoretical approaches have been performed to an- 
alyze ferromagnetism in DMS, this includes mean field 
theory_f|-|5| , spin wave theory || , first principle calcula- 
tions and Monte-Carlo simulations ]Tcj] - In contrast 
to most of the theoretical work, we present a theory is 
able totreat disorder in a more realistic manner (beyond 
coarse graining). Our theory includes quantum fluctua- 
tions within RPA and disorder is treated within CPA. It 
should be stressed that in our approach the spin impuri- 
ties are treated quantum mechanically. 

We start with the following minimal Hamiltonian, 



H ~ tij C ia C j<? 



^ JiSi 



(1) 



The first term stands for the tight binding part of the 
itinerant free electron gas, £y = t if i and j are nearest 
neighbor or otherwise. The second term is the exchange 
between localized impurities spin and itinerant electron 
gas, Ji are random variables: J; = J if site i is occupied 
by a Mn 2+ ion or 0. The operator Si — (y itx {l/2a a p)Cif) 
is the spin operator at i of the itinerant electron gas and 
Si is the spin of the magnetic impurity. 
Let us define the Green's function, 

Gt(t) = -ie(t)([S+(t),Sr(Q)}) =« S+;S- » (2) 



We write the equation of motion and use Tyablicov 
decoupling |ll[] (equivalent to RPA) which is suitable for 
ferromagnetic systems. It consists in closing the system 
by approximating the higher order Green's function <C 
S-sf; Sj (S*) < s+; Sj >. In this approximation, 
we obtain in frequency space reads, 



(u + J i (a'))G%-(u,) = 2(S;)6 ij + 



(3) 



(s z ) is the magnetization of the itinerant electron gas 
and (Sf) the magnetization of a magnetic ion at site 
i. It is convenient to rewrite the new Green's function 
which appears in the right part of the equality in the fol- 
lowing form, < s+;Sj »= £ £ fcg e ^T h j +q ' k , where 



T k+ q ,k =< C t + ^ cu .£- > . We obtain, 



-\k-\-q,k 



f(k,q,u)J2ljie- iqRl G 



where, 



f{k,q,w) 



(( n k+q,t) - ( n k,l)) 

u — (eh — £k+ q ) + cJ (S\) 



(4) 



(5) 



( n k,a) is the occupation number of (k, a) state, c is the 
impurity concentration, (S\) is the averaged magnetiza- 
tion of Mn 2+ , and denotes the hole's dispersion. In- 
serting both eq. (f|) and (||) into eq. (||) we immediately 
find, 



i 



where the T dependent locator gi is defined as, 

2(5?) 



LO + Ji(s z ) 



(6) 



(7) 



4>u = jJiJiXai^) an d X?i( w ) is the Fourier transform of 
the polarized susceptibility x°(g',a;): 
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(6A 



(8) 



Note that eq. siiZZ contains the disorder through 
4>u and It is also interesting to mention that the 
previous equation can be interpreted as the propagator 
of a free particle moving on a disordered medium, g$ is the 
random on-site potential and <pu the long range-hopping 
terms. Note also that 4>u is energy dependent through 
X° ; (w). To solve the problem we have to calculate in a 
self-consistent manner (s z ) and (rik M ) which appear in 
eq. (||). For that purpose we have to write the equation 
of motion for the Green's function Kij l<7 =<C Ci l<T ; ct a 3>. 
After decoupling we get, 

{u-^z a J i (Sf))K ii , a = t 



(9) 



One can recognize the propagator of the Anderson 
model, with on-site random potential depending on the 
spin a: — \z a Ji(Sf). Since in our model the po- 
tential is temperature dependent through (Sf) then suf- 
ficiently close to T c we will always be in the metallic 
regime kfl e 3> 1 l e ~ rjr^yp ■ This is in contrast 

with the standard Anderson model where the impurities 
are static. Equations @ and (g) (a = ±1) provide a 
closed system of equations which have to be solved self- 
consistently within CPA. 

The simplest is to start with eq. (||). Indeed, it is 
straightforward to get the solution with the standard 
CPA since it contains only diagonal disorder. The av- 
eraged Green's function is 



1 



(10) 



uj - e(k) - T, a (uj) 

where the self-energy is, 

E„(u) = V a - [e A ,„ - V a {u))K™{u){e B ^ - (11) 

where ca,ct = \ z oJ (S\) ■, ts.a — and V a is the average 
value V a =\z<rMS%) and K°° = ± £ g .- e(g) - W 
The self-energy E CT (a;) can re-expressed, 



= V a [l + (1 - c)-z a J{S A )K™{u)] 



(12) 



We see that when T -> T c , E CT (u;) -> £^ CA M = V a . 
Thus, in the framework of our decoupling scheme, close 
enough to T c the CPA for eq. (10) reduces to the VCA. 

The final step of the calculation consists in solving 
eq. (](]). In order to provide analytical form for T c we 
use the similar approximation (VCA) for G~l~ as done 
above for 2£y i<T . We expect this approximation to be rea- 
sonable in the limit of both low impurity concentration 
and low density of itinerant carriers. To get the aver- 
aged Green's function, we use the well-known Blackman- 
Estcr ling-Beck formalism Jl3|,[l4[ . By contrast with stan- 
dard CPA, this approach is suitable for non diagonal dis- 
order problems. It is based on a 2 x 2 matrix Green's 



function formalism for binary alloys using locator expan- 
sion. Within VCA approximation, one gets for the aver- 
aged Green's function of an atom of type A, 



G+- 



9a -ca(k,uj) 



(13) 



where a(k,uj) = |j 2 x°(fc,w). 

Note also that since the Ga atoms have no magnetic 
moment, it implies G~^~(k,oj) = gs = 0. The Mn 2+ ion 



propagator can be rewritten, G A (k, ui) 



E = 



E-E(q) 



The dispersion E(q) is solution of, 



E(q) = -J^ + ±J 2 X (q,E(q)(S z )) 



, where 



(14) 



According to ref. |TJj the magnetization can be ex- 
pressed in the following form, 



(Sl) = 



\2S+1 



(5 + 1 



2S+1 



(15) 



whCTe <f>= l£q e^)-r 

When T -»■ T c , </> = Mki 
T c the standard RPA form 



E q eW)- Thi s implies for 



rp 1 S(S 



(16) 
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This expression is similar to the one obtained in the 
clean limit for the Kondo Lattice Model Jl^,|l^]. In the 
vicinity of T c the dispersion E(q) is, 



1 J 2 I 



E ^ = ^T2^--q^-V l < 



2k f 



2k, 



)) (17) 



Note that, below T c , the eq. ( |14[ ) should be solved nu- 
merically in order to get E(q) as a function of the tem- 
perature. This is required to calculate (5 Z ) and (s z ) as 
function of T. According to eq. (|l6|) T c is given by, 



5(5+1) J 2 c. 1 



24tt 2 



t K N ^ 



C(q,k f ) 



(18) 



where we define C(q,k f ) = \{k f - \\k) - £]I„(g*£)). 

This implies that T c is proportional to J 2 and to the 
effective mass (1/i). The dependence on the hole con- 
centration is only contained in C(q,kf). We define the 
hole concentration as rih — 7 c where 7 < 1. This is the 
simplest way to take into account the presence of As- 
antisites p8[ . In Fig. 1 we show the variation of T c as a 
function of 7. 



2 



250 



200 



VCA 

Stiffness approx 
Mean Field 




- 5 log(T)° 




FIG. 1. T c as a function of 7 for c = 0.05 and ^- = 10.5 eV. 
The continuous line represent the VCA calculation, the 
dashed line corresponds to T c within the mean-field approxi- 
mation, and the dotted line is obtained by approximating the 
dispersion by E(q) = Dq 2 , D is the spin stiffness. The inset 
shows log(T c ) versus log (7). 

We observe that in the low hole concentration regime, 
T c agrees very well with the mean field result (this 
is more clear in the inset log-log plot). In the mean 
field regime the magnon excitation spectrum is disper- 



lim„ 



sionless: E MB (q) — mi, r 
kf = (37r 2 7c) 1 / 3 . In this limit, 



1 , 3 



E(q) 



J 



1 

Sw 5 



r« = 5 (/i(s + i)r» 



1/3-4/3 



(19) 



When increasing 7, T c strongly deviates from the mean 
field results and shows a broad maximum. Such a max- 
imum was also observed in ref. 0. By further increase 
of 7 the Curie temperature starts to decrease As 
we observe it from Fig. 1, for very large 7 T c agrees very 
well with the case where the magnon spectrum is approx- 
imated by E(q) = E stlf t (q) = Dq 2 where the stiffness D 
is given by D = j^-? ^- -rj , this regime is denoted "stiff- 



48tt t t kf 

ness" regime. In this regime we find, 



1 



1 



144(18tt 4 ) 1 /3 
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l)Z! 7 -l/3 c 2/3 



(20) 



The existence of a maximum can be understood in the 
following way: Like in the RKKY situation [pofl , the ex- 
change oscillate with typically length scale l osc oc 
Thus it is expected that when the length scale gets suffi- 
ciently large (larger than the average distance between 
impurities) some Mn-Mn bonds are coupled antiferro- 
magnetically. The induced frustration has for immediate 
consequence a decrease of T c . In Fig. 2 we illustrate the 
previous discussion by showing the dispersion as a func- 
tion of k/k c where k c is chosen in order to conserve the 
volume of the Brillouin zone (v = (27r) 3 ). The results are 
shown for the 3 different regions: "mean field" , "inter- 
mediate" and "stiffness" regime. We observe that in all 



cases the dispersion goes to (when q — » 0), as expected 
when the Goldstone theorem is fulfilled. 
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FIG. 2. Magnon dispersion for c = 0.05 in the 3 different 
regime, "mean field" (continuous line), "intermediate" (dot- 
ted line) and "stiffness" regime (dashed line). E(q) is rescaled 
by a factor A = ^(3tt 2 7) 1/3 

In Fig. 3, we show the region for which T c reaches 
its maximum as a function of c (T™ ax (c)) and the re- 
gion where MF formula provides a good approximation 

rpMF 

— I < 0.1. First we see 



-kf where for T c , it corresponds to 



that the region of validity of the MF result (dashed area) 
corresponds to a very narrow region typically 7 < 0.05. 
A good approximated value of the 7 for which T c is max- 
imum can be obtained by taking the intersection point 
between the MF and "stiffness" values. This leads to 
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FIG. 3. The dashed area represents region where 
mean-field result for T c is valid, the symbols are calculated 
points and the dashed line a fit. The continuous curve rep- 
resents values of (7,0) for which T c is maximum (T™ ax ) . In 
the inset we have plotted T™ ax as a function of c assuming 
^ = 10.5 eV 

So far, for our discussion we did not have to specify 
the values of the parameters t and J. In order to check 



3 



the validity of our theory we compare our results with 
available experimental data. GaAs is known to have a 
fee structure with a lattice constant ao ~ 5.6A For sim- 
plification in our calculation we have assumed a simple 
cubic structure thus the lattice constant which has to be 
taken in our calculation is a\ = ^773 in order to conserve 
the volume for the unit cell. By also assuming an effec- 
tive mass for the holes m = 0.5 m e one gets t = 0.63 eV. 
The remaining free parameter J will be chosen in order to 
fit the experimental data of ref . || . For that purpose we 
calculate 7 for each sample according to the measured ex- 
perimental values of the hole concentration given in Fig. 
2 of ref. The results are depicted in Fig. 4. As it can 
be seen we find a very good agreement with the experi- 
mental data if ^- = 10.5 eV, this implies |J| = 2.58 eV 
p2| . Note that the deviations observed at low c are due 
to the uncertainty on the hole concentration value (see 
the huge error bars in Fig. 2 of ref. 0). 
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FIG. 4. T c in Kelvin as function of c. The full square 
corresponds to experimental values taken from ref. Q|. The 
full circles represent the calculated values, the 7's where also 
taken from the same reference. 

From the experimental measurements, there is no clear 
consensus concerning the correct value of this parame- 
ter. Indeed, recent core level photoemission has provided 
J = -1.2 ±0.2 eV Q. Whilst, from Magneto-transport 
measurements a value of | J| = 2.4±0.9 eV was suggested. 
P, pl|. And within first principle calculations Sanvito et 
al. have found J » -4.65 ± 0.25 eV. In order to 
proceed to a better estimation of the parameters J one 
should compare theoretical calculations with other data, 
for instance transport measurements data [ p4| . However, 
it is interesting to note that the band splitting at T = OK 
(A = JcS) obtained within our calculations agrees with 
the experimental value reasonably well fl25| . In the inset 

t2 

of Fig. 3, assuming =r- = 10.5 eV we show T c as func- 
tion of c taking 7 on the line of "maximum of T c " . For 
instance if c w 0.1 and 7 w 0.2 a T c of order 230 K can 
be reached. 
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FIG. 5. Normalized magnetization as function of T/T c . 
The experimental data are taken from [B]. 

Let us proceed further on and compare the calculated 
magnetization with the measured one. In the experimen- 
tal data the concentration of Mn 2+ is c = 0.053% and 
the parameter 7 is estimated to be 0.3 (see ref. [^lj) In 
Fig. 5, we show the magnetization as function of the T 
where m* = 0.5 m e and J = 2.58 eV. . We observe that 
for sufficiently high temperature T < 0.5 T c there is a 
very good agreement with the measured Mn 2+ magneti- 
zation. When decreasing T some deviation appears, this 
suggests that the VCA treatment is not good enough in 
this region, which was expected. 

To conclude, we have presented a general theory for 
carrier induced ferromagnetism in DMS. Our approach 
allows to treat the disorder beyond simple coarse grain- 
ing within full CPA treatment. It goes beyond mean- field 
and includes quantum fluctuations in the RPA approx- 
imation. We have shown that, within our decoupling 
scheme and sufficiently close to T c the CPA for the itin- 
erant gas reduces to VCA, which allowed us to provide 
analytical results for T c in the low impurity concentration 
and hole density regime. We have also discussed its de- 
pendence on the hole concentration. We have also shown 
that the mean field approximation is only valid for very 
low carrier concentration. Additionally, for illustration 
of our theory a comparison with available experimental 
data on Ga\- c Mn c As was done. We find a very good 
agreement with the experimental results assuming a sin- 
gle band for itinerant carriers and a large exchange con- 
stant J = —2.58 eV. Finally, this work provide a good 
starting point for higher decoupling scheme. 

Note added: After this work was completed we be- 
came aware of Yang et al. comment fl7|| . By analogy 
with the Kondo Lattice Model |l6| (no disorder) they 
proposed a similar expression for T c to the one derived 
in eq.(|l6|). 
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